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Abstract 

We studied meson-meson interactions using the 't Hooft model, which represents QCD in 1 + 1 
dimensions and assumes a large number of colors (N c ). The dominant interactions in this large 
N c limit are generated by quark exchange. Our results show that QCD in 1 + 1 dimensions allows 
the realization of a constituent-type quark model for the mesons, and generates a scalar "<r"-like 
meson-meson resonance, whose effective coupling and mass are determined by the underlying QCD 
dynamics. These results suggest an interpretation of the lightest scalar mesons as qqqq systems. 

PACS numbers: 13.75.Lb, 13.75.-n, 12.38.-t 
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I. INTRODUCTION 



At low energies, systems which are due to, or interact through, the strong nuclear interaction 
may be described by effective field theories, examples of which are Quantum Hadrodynamics 
and Chiral Perturbation Theory [2]), or simple constituent quark models [3|. With 
more or less phenomenological content, those frameworks are in general very successful, but 
there are still questions not completely answered. 

For instance, whether the nature of the broad a meson corresponds truly to a simple 
quark structure, or to a resonance in the meson-meson dynamics, or to an unusual quark 
structure as a meson-glueball combination, or even to some combination of these, is still 
an open issue. Models of QCD did not yet resolve this question quantitatively. In general 
for the light scalar mesons, arguments have been advanced for the importance of a qqqq 
component |4| at short distances, compatible with a dominant meson-meson component 
at large distances. A very recent lattice calculation [5( also indicates that a qqqq bound 
state may exist, just below threshold in the non-exotic channel of pseudoscalar-pseudoscalar 
s-wave scattering. 

In this work we apply the 't Hooft model jfj, a formulation of QCD in 1 + 1 dimensions 
and in the large N c limit, to the meson-meson scattering process. The 't Hooft model has no 
physical gluons, thus it includes no glueballs. Furthermore, due to the large N c limit, quark 
exchange dominates gluon exchange and only a finite number of diagrams contribute. Since 
the finite sum of regular contributing diagrams is unable to produce a pole at real energies, 
no meson-meson bound states can be produced. 

However, while meson-meson bound states are excluded in this framework, complex en- 
ergy resonant states should still be possible. It is the purpose of this paper to calculate 
the meson-meson scattering amplitude based on meson-gg vertex functions obtained in the 
't Hooft model, and to look for low-lying resonances. As we will show, the 't Hooft model 
indeed supports the existence of meson-meson resonances, suggesting the relevance of the 
qqqq structure for the light scalar mesons. 

A calculation of 7r — 7r forward and backward scattering in the Dyson- Schwinger, Bethe- 
Salpeter approach and in the rainbow-ladder approximation was presented in Ref. |7J . It uses 
an effective qq interaction and incorporates features of QCD. Based on this approach, a scalar 
meson emerges as a resonance in it — ir scattering. However, the calculation is performed 
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using the Euclidian metric. In our work, while simplifying the problem by working in 1+1 
dimensions, we use the Minkowski metric througout. 

Clearly, a model in 1+1 dimensions is limited in its scope, and one has to be very cautious 
when comparing its results to phenomena in the real world. Nevertheless, for kinematic 
conditions of scattering processes with a negligible component of the momentum transfer in 
the transverse direction, we may conjecture that the 't Hooft model, and thus the calculation 
presented here, has the main features of realistic microscopic QCD, and that its results are 
valid at least qualitatively. 

Section |H] reviews the model and introduces the calculational framework. Section IHII 
presents the input for the description of the meson-meson scattering transition amplitude. 
Section HVI shows the results, and section IVI presents the conclusions. 

II. FORMALISM 

A. The 't Hooft model and the choice of gauge 

This work models the meson-meson interaction using the 't Hooft model. Here we review 
briefly the dynamics of this model by starting with the corresponding Lagrangian. Subse- 
quently, we write the equations that we solved for the one body (quark propagator) and two 
body (quark-antiquark bound state) problems. 

The QCD Lagrangian is 

C = - l -Tr [G^G^] + q [iDyrf -m )q, (1) 

with the notation 

a 2 > 

Gpv = d^A v - d u A^ + ia [A^ A v \ , 

D^= 3^ + iaA^ (2) 

where A^ are the gluon fields with the Lorentz index /i and the color index a, the A a 's are 
the generators of the SU(N) color group, G^ u is the field tensor, q is the quark field, m is 
the bare quark mass and a is the quark-gluon coupling strength. Following 't Hooft p , the 
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coupling strength a depends on the number of colors in the following way: 

Qt=-p=. (3) 
Introducing for an arbitrary two- vector b the light cone variables: 

6+ = 7f (6 ° + 6l) ' 

b- = -^(b - b 1 ) , (4) 

the scalar product of any two vectors a and b becomes a + 6_ + a_6 + , and the derivatives 
correspond to 

o d i ! i <9 a x 

9 - = 71— = - d ) = "^(^3T + 



<9s + a/2 a/2 5a; cte 1 

a+ = £ = >° + *> = 7i<ip-£r>- < 5) 

In the same way, the + and — components of the 7 matrices are defined. The anticommu- 
tation relations read 

{7-, 7-} = {7+, 7+} = 0, 

{ 7+ , 7 _} = 2. (6) 

In the light cone variables, the nonvanishing components of the field strength tensor are 

G+_ = = <9+A_ - d„A + + ia[A_ } A+] , (7) 

so that the free gauge field Lagrangian becomes 

We choose to work in the light cone gauge, A_ = 0, so that the commutator contained in 

the field tensor disappears. Moreover, since we consider only two dimensions, we do not 

have any physical gluonic degrees of freedom. In addition, because of the gauge condition, 
there is only one degree of freedom left. Consequently, the gluonic field is not a dynamical 
variable and does not couple to ghosts any longer. 

In this parametrization the Lagrangian (0) becomes 

1 



£=-IY 



(cL Af ) + q (i<9+7_ + id^ + - a^-A + - m ) q . (9) 



Before quantizing the theory given by this Lagrangian, we calculate the gluonic field. The 
equation of motion related to this field is 

(^-) 2 A + = -«g 7 _g. (10) 

The solution of Eq. ([TUj) is 

= -<./*^,«-h-^.«-)Cto + -x + ), (ID 
where the Green's function Q is given by 

G(y+ ~ x+) = \y+ ~ x+\ + ci(y+ - x+) + c 2 . (12) 

The coefficients C\ and c 2 are free parameters. This means that the gauge condition did not 
eliminate all the redundant degrees of freedom, just as the Coulomb or Lorentz gauge do not 
determine uniquely the photon propagator in QED (Gribov ambiguity). We can therefore 
set the coefficients c\ and c 2 equal to zero in order to simplify our calculations. 

The Fourier transform of the Green's function ()12|) gives the gluon "propagator" , or more 
precisely the momentum dependence of the effective quark-quark interaction: 

D(k) = £>(*_) = i - 6{k-)P jT ^ . (13) 

The second term in Eq. ()13j) was first considered by Gross and Milana [8| , in the different 
context of a quasi-potential two-body equation for the quark-antiquark system. It makes 
the potential A + finite everywhere. 

From this point we proceed to solve the one body equation for the quark propagator. 



B. Quark Dyson-Schwinger equation 

The (undressed) fermion propagator is 

So{k)= 2M--mg + i6 ' (M) 
and the quark-gluon interaction is 

— iV = —ia^- . (15) 
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FIG. 1: The quark Dyson-Schwinger equation. The curly line represents the strong interaction, 
the thin line the unperturbed quark propagator, and the solid line represents the dressed quark 
propagator. 

We determine the dressed single quark propagator S(p) using the (one body) Dyson- 
Schwinger equation (DSE), 



which we show also graphically in Fig. ^ 

Since for every internal loop there is a factor of a 2 = g 2 /N c , and a multiplicative factor 
of N c , the color dependence disappears. The vertex corrections and the quark-gluon vertices 
do not have a multiplicative factor, being supressed in the large N c limit. Therefore, in 
this limit the rainbow approximation (undressed vertices and the absence of the quark loops 
from the gluon propagator) is justified Q]. 

In Eq. (fTHjl d 2 k = dk-dk + , and since D{k) does not depend on k + , the principal part 
V J d 2 k does not depend on p + . This allows the following parametrization of the full quark 
propagator, 



S(p) = S (p)+ig 2 S(p)[V J 



d 2 k 



D(k-p)^S(k)^]S ip), 



(16) 



(2tt) 2 




(17) 



where the self-energy contribution E, originated by V J d 2 k, depends only on k_: 




2k_k + - £;_£(£;_) - ml + ie ' 



(18) 



Performing the k + integral first one obtains 




(19) 



Substituting this result back into Eq. (fTSj) one finds that 




(20) 
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FIG. 2: The Bethe-Salpeter equation for the qq bound state. 



Using Eq. (jlHjl for D(k- — p~) and performing the integral we get 



E(p. 



7T£>_ 



This, in combination with Eq. (|T7|) results in 

(P+ + ^f)T- + + ™o 



5(p) 



2p+P- - (mg - ^ - 



Note that the masspole has been shifted: 



^ - ...2 _ 

7T 



m — > m = m — 



(21) 



(22) 



(23) 



Having obtained the dressed propagator, we are ready to proceed to the next stage, namely 
the calculation of the q-q bound state. 



C. Two-body bound states 

In the following, we label the dressed quark mass by m\. As for the antiquark (which 
might have a different flavor), we label its dressed mass by m 2 . The total momentum of the 
bound state is denoted by r and the momentum of the quark by p. The momentum of the 
antiquark is then r — p. 

The bound state wave function T(p, r) is given by the Bethe-Salpeter equation (also shown 
graphically in Fig. EJ), 



f d 2 k 

T(p,r) = ig 2 V / j—^D(k-)^-S 2 {p + k - r)T(p + k,r)S 1 (p + k)j- 



(27T) 



(24) 



where S\ and S 2 are the quark- and the antiquark propagators, respectively. With the 
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substitution T{p,r) = 7_?/>(p, r) p Eq. (|2*3)l becomes 
tfj{p,r) = i(2g) 2 V 



2 I d 2 k D(k )(p + k)_(p + k — r)_ip(p + k,r) 



(2tt) 2 [{p+k) 2 -m 2 ][{p+k-r) 2 -m 2 



(25) 



The equal x_ wave function is defined in the following fashion: 



(p{p_, r J)= dp + — 2TT7 V2 2T 26 



By substituting this into Eq. (J2o"|) one gets 
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^(p, r ) = -^-P / dk-D(k-)ip{p- + fc_, r_) . (27) 

Note that ip(p,r) does not depend on p + . Multiplying both sides of the former equation by 
\p-(p — r)J\/[(p 2 — m\)((p — r) 2 — m 2 ,)] and integrating over p + we consider the two poles 

2 2 

in the complex p + -plane, namely k\ = ^ — iesgn(p_) and h% = 2( - p m _! r -\ — ^esgn(p_ — r_). 
If both of them are in the same half plane the integral over p + is zero, because the sum 
of the two residues is zero. If the first pole is in the upper half-plane and the second one 
is in the lower half-plane (which means that p_ < and p_ — r_ > 0) the integral is 
2m/ (ki — k 2 ) = 2m9(—p-)9{p- — r_)/ (ki — k 2 ). If the second pole is in the upper half-plane 
and the first one is in the lower half plane the integral is —27ri9(p_)9(—p_ + r_)/(/ci — k 2 )- 
Combining these two cases one obtains 

r) = f { - pJ) %- r - ) -J"' P - )B{r - pJ) ^ r) . (28) 

2p_ 2(p_-r_) r + 

Whenever the combination of the 9 functions does not vanish, it is easy to invert this relation: 

ip{p, r) = — (8(-p-)6(p- - r_) - 0(p_)0(r_ - p_)) — — - — - r + cp(p, r) . 

7U \2p~ 2(p^ — r_) / 

(29) 

Whenever this condition does not stand, we have to use Eq. (J27)) to compute ip from Note 
that ip has been normalized to 1/ \fW c in order to get the correct charge normalization. 

In order to determine ip, we transform Eq. f)26j) into a form suitable for a numerical 
calculation. The 9 functions limit the range of p_ to < p_ < r_, and for real particles 
only positive values for r + must be considered. After some algebraic manipulations, Eq. 
(J2HJ) becomes 

^^e +T ^),(-)-p/;^^|^, (so, 

where the following notation was introduced: 



2 27rr + r_ 7rm? 7rmn P- L . . 

V = - 2 — , «i = — J-, a 2 = ^^, x= — , y= — . (31) 
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FIG. 3: The bare mass (solid line) and the dressed mass (dashed line) of the first quark, as a 
function of the strong coupling constant g, with the constraints that the pion mass is 140 MeV 
and the ratio moi/m-02 = 3/4. 

We solve the integral equation ()30j) numerically. The wave function is expanded in cubic 
splines (since the wave function ip is defined only in the range between x = and x = 1, the 
boundary condition that they vanish at the limits of this interval is imposed). The resulting 
linear matrix equation for the expansion coefficients was solved with standard eigenvalue 
routines. 

In the limit mm = m 02 = 0, Eq. (fHUj) yields a ground state of zero mass, and thus is 
consistent with chiral symmetry. To generate a solution that is related to the pion in the 
real world, we searched for a bound state solution of Eq. (J3*U|) with a mass of 140 MeV. To 
obtain such a solution, we varied the bare mass moi of one of the quarks. For simplicity, 
the second quark mass was not taken as an independent free model parameter but 
determined by assuming a fixed mass ratio =3/4, which lies within the accepted range 



between 0.2 and 0.8 (lOj]. The coupling parameter g and the dressed masses m± and m 2 were 
adjusted accordingly, through Eqs.tjHOJ), (EIQ) and (|2l3|) . 

We represent in Fig.|3]the values of the bare quark mass m i, as a function of the coupling 
strength g, which generate a bound state with a mass of 140 MeV. The bare masses are found 
to depend linearly on g. The slope and the y-axis intercept of the numerical straight line on 
Fig. El are easily determined through a fit, with the result 




m i = 60.57- 0. 178 # MeV. 



(32) 
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With the help of Eq. (J2*3|) we can also predict the dependence of the dressed masses on g 
from the curves for the bare masses. Therefore, in the 't Hooft model the dressed masses 
are given as the following functions of g: 

a 2 

ml = (60.57- 0.178 gf - — MeV 2 

7T 

a 2 

m\ = (80.76 - 0.24 g) 2 - — MeV 2 . (33) 

7T 

We can also determine the largest value of g, such that each dressed mass is physical, i.e., 
not imaginary. For the first flavor this happens at g — 81.64 MeV, while for the second 
flavor at g — 100.8 MeV. The first value is therefore the largest possible for g, such the 
't Hooft model may support a constituent quark model interpretation, where the dressed 
masses correspond to constituent quark masses. 



III. MESON-MESON SCATTERING 



In this section we consider the meson-meson elastic scattering amplitude. We continue 
to assume two different flavors for the quarks, whose dressed masses are m\ and m 2 , and we 
consider the lowest qq bound state only. 

The diagrams that dominate in the large N c limit are the quark exchange box diagram, 
represented in Fig. El and the quark exchange crossbox diagram, represented in Fig. 0D In the 
center of mass system, the momenta of the ingoing mesons are P = (P°, P 1 ) = (y/ fi 2 + p 2 ,p) 
and P = (P^P 1 ) = (a/ fi 2 + p 2 , —p), where /i is the mass of the meson and p the relative 
momentum. The outgoing particles then have the same (but interchanged) momenta. 

Both diagrams are symmetrized in terms of the outgoing states. Similar diagrams which 
are obtained from the former ones by interchanging m\ and m^ in the intermediate state are 
also considered. There are a total of eight diagrams which were calculated. When there is 
only one quark flavor, one does not need to interchange the two masses and there are only 
four diagrams. The sum of these diagrams is proportional to 1/N C . 

As for the gluon exchange diagrams, such as in Fig. El they are suppressed in the large 
N c limit by a factor of 1/N C compared to the quark exchange terms. 

Since the vertex function is independent of the + component of the relative momentum, 
the momentum integral in the loops of Figs. |U and El is simplified: we can first integrate 
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FIG. 4: The quark-exchange box diagram. The vertex functions are represented by filled circles, 
and the mesons by double lines. The power counting from the vertices is explicitely shown. An 
extra factor of N c comes from the color summation in the internal loop. 
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FIG. 5: The quark-exchange crossed box diagram. As in Fig. the power counting factors are 
explicitely shown. An extra factor N c comes from the color summation in the internal loop. 

the propagator products over k + analytically, and then perform numerically the second 
integration over which includes now the vertex functions. 

As an illustrative example, we demonstrate the calculation of the box diagram (Fig. @J 
in greater detail. 

The corresponding scattering amplitude is 

M hox = / dk_ / dk + ij(-k,-P)ij(k,P)^(P + P-k,P)iP(k-P-P,-P) 

JOD JOO 

1 1 

^ 2 ~ 2 ~ 

k + ~ St + ie s S n ( k -) k+- P-- m *p + ie sgn(A;_ - P_) 
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FIG. 6: The gluon exchange diagram. Power counting is shown as previously. For this purpose 
the gluon line is represented as two parallel quark lines. An extra factor N c comes from the color 
summation in the internal loop. 
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k + - P + - P + - 



2(fc_-P_-P_) 
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+ ie sgn(/c_ — P_ — P) 
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The propagators have four poles: 



+ ie sgn(/c_ — P_) 



(34) 



fci 
ki 

k 3 



m 



.^ k L -iesgn(k^), 



2(A;_ - P_) 



iesgn(A;_ — P_ 



P+ + P+ + 



ml 



2(k_ 



P-) 



— ie sgn(A;_ — P_ — P_) , 



h = P+ + 



iesgn(/c_ — P_) . 



(35) 



2(fc_ - P_) 

In order to perform the k + integration, one needs to close the contour in the complex plane 
and consider the residues of all poles inside the contour. There are 16 different possible 
combinations of signs of the imaginary parts of the poles. Some of these cases can be 
excluded, because they correspond to values of fc_ which make the k + integral vanish. 

For instance, a pole k\ in the upper half plane implies that the pole A; 2 cannot be in 
the lower half-plane, otherwise one would have A;_ > P_ = (y/ fi 2 + p 2 + p)/ V2 > 0, in 
contradiction with the initial hypothesis k- < 0. Likewise the poles k s and k 4 cannot be 
in the lower half-plane either. Therefore, if k\ is in the upper half-plane, the other 3 poles 
are also in the upper half plane. This would imply the k + integral to vanish, since one may 
close the contour below the k + axis. Therefore we can exclude the case when k\ is in the 
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upper half plane. 

After a detailed analysis one finds that there are only three cases that have a non- 
vanishing contribution to the integral: (i) only k 3 is in the upper half plane, (ii) the poles 
k 2 and k 3 are in the upper half plane, (iii) only k\ is in the lower half plane. As for case 
(i), it implies k_ > P_ and < P_ + P_. Under these circumstances, ip(P + P — k, P) = 
^(k 3 -k 2 )\ e=0 <p((P+P-k)^/P^ a nd^(k-P-P,-P) = %(k 2 -k 1 )\ e=0 <p({P+P-k)-/P-), 
due to Eq. while the other vertex functions have to be evaluated using Eq. (|27jl. The 

contribution from case (i) becomes 

A^boxli = [ - ) /. dk 



it J Jp_ (h-k 4 )\ e=0 
x y( (P + P p ~ k) - M (P + p~ k) ~ M-k, -PMk, P) ■ (36) 

This integral is computed numerically. We treat the other two cases in a similar fashion. 

We note that the amplitude of Eq. (j3T?j) vanishes in the chiral limit, which supports that 
the ground state solution for the qq system has features of the pion. Indeed, since in the 
chiral limit Eq. (J30|) gives \i = 0, one has P_ = ^(Va* 2 + P 2 — p) = 0, implying that the 
two integration limits in (pffij) coincide, and therefore the amplitude vanishes. The same can 
be shown for the other terms of the amplitude not explicitly written here. 

It is worth mentioning that we implemented stability tests of the numerical results against 
the number of gridpoints, the number of splines and the singularity regulator e (the k + 
integral of the propagators is singular). These checks proceeded by imposing the following 
criteria: doubling each of the mentioned parameters, the relative change in the results 
should be less than 1% . Convergence is typically obtained for 440 gridpoints, 40 splines 
and e = lO" 2 . 



IV. RESULTS 

The numerical calculation of the meson-meson scattering amplitude (section IIII|) starts 
with the evaluation of the two-body quark-antiquark wavefunction from the Bethe-Salpeter 
equation (section III CJ) . In turn, this demands as input the bare quark masses and the 
quark-gluon coupling constant g (section IlIBjl . 

We constructed four representative models which correspond to four different choices 
of the quark-gluon coupling constant and bare quark masses. They are subjected to the 
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Model I 


Model II 


Model III 


Model IV 


g 


1 


20.1 


80 


2500 


moi 


60.0 


57.0 


46.4 


6.0 


m 02 


80.0 


76.1 


61.8 


5.0 


mi 


60.0 


55.9 


10.6 




m 2 


80.0 


75.2 


42.2 


- 


m r (= E R ) 


280.0 


282.4 


280.7 


300.0 


m 


4.4 


68.6 


80.7 


139.3 


r 


0.07 


16.4 


23.2 


64.7 



TABLE I: Constants and resonance parameters for models I— IV. The first five lines show the quark- 
gluon coupling constant, the bare the dressed quark masses (the latter are unphysical in model IV 
and thus not included). The following lines are the results of the s-channel resonance fit using Eq. 
(|37j). All parameters are in MeV. 

constraint that the bound state mass of the qq system is the pion mass 777,^=140 MeV. 

The parameters denning the four models are shown in Table HJ The organization principle 
for these models is that, going from model I to model IV, the quark-gluon coupling constant 
increases and the quark masses decrease. 

In models I and II, the sum of the dressed quark masses is close to the real pion mass. 
Consequently, they can be interpreted as constituent quark masses in a constituent quark 
model for the "pion", generated by QCD within the 't Hooft model assumptions. However, 
this correspondence breaks down for large couplings, as seen for the model considered next. 

In model III, the value of the coupling g is chosen slightly below the maximum value 
determined in section 2 for a constituent quark model interpretation, while in model IV that 
maximum value is exceeded substantially. In this last case, the dressed masses are imaginary 
and thus not physical. 

Not surprisingly, the bare quark masses differ from the up and down current quark masses 
of QCD for all models. This is a known artifact of the 't Hooft model: from Eq. (|2*3"|) one can 
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FIG. 7: Bethe-Salpeter qq wave functions for models I-IV. x = ^ is the momentum fraction (or 
Bjorken variable) defined in Eq. (|31|). 

see that the mass shift due to the dressing decreases the quark mass, instead of increasing 
it as in 3+1 dimensional QCD. 

Figure [7| shows the obtained Bethe-Salpeter qq wave functions for each model. Clearly, the 
description of the quark-gluon vertices varies considerably. The wave functions are strongly 
peaked around x close to 0.4 for small values of g, while for larger g they become broader 
and more and more constant. This behavior can be easily understood: larger values of g 
cause stronger attraction between quark and antiquark, leading to a tighter bound state and 
therefore a more spread out wave function in momentum space. 

As described in the previous section, for each of these models we calculated the meson- 
meson scattering amplitude, the squares of which are displayed in Fig. |H1 Since we are 
primarily interested in their structure, the amplitudes have been scaled such that the maxi- 
mum of their absolute squares are equal to 1. In all four cases, we find a resonance structure 
close to threshold. 

This feature could be a sign for the existence of a q 2 q 2 bound state, for which the lattice 
calculations of Ref. jjj found indications. We remind the reader that working in perturbation 
theory we cannot generate a bound state directly, but it would be interesting to see if in a 
non-perturbative extension of our calculation, such a bound state would also emerge from 
the 't Hooft model. 
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FIG. 8: Absolute squares of meson-meson scattering amplitudes for models I - IV close to threshold 
energies. A resonance structure is clearly visible. The solid lines are obtained using the 't Hooft 
model, the dashed lines represent fits to a simple resonance model. The amplitudes are scaled to 
1 at their maximum. 



On the other hand, the experimentally observed resonances have energies well above 
threshold. It may be necessary to include gluon exchange (not considered in our calculation) 
and higher order quark exchange terms in the expansion in powers of l/iV c in order to achieve 
a description resembling more closely the real world. 

In order to determine approximately the position and width of the resonance, we compare 
the 't Hooft model amplitudes to a simple resonance model. We calculate the amplitude for 
an intermediate s-channel resonance at tree level, the absolute square of which is 



\M? = S* ,_J. 2+ . , (37) 



where vr? r and mf are the real part and the imaginary part of the square of the resonance 
mass, respectively, and g is the effective meson- meson- "cx" coupling constant. We then fit 
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the parameters of this simple model to best reproduce the 't Hooft model results. 
In the non-relativistic limit, (|37j) becomes 

1 



\M\^t~ ^— ^r- (38) 



We can compare (Jag) with the well-known Breit-Wigner form, 



i*i'"^5W < 39) 

where ^ is the resonance energy, T is its width, and g the corresponding coupling strength, 
and simply read off the relations T 2 = mj/rri 2 . and Er = m r . 

The parameters obtained in this way for all models are also displayed in Table HJ We 
verified in an independent fit directly to the Breit-Wigner form (|39|) that the parameters are 
not significantly altered in the non-relativistic limit. For comparison, the mass of the "real" 
sigma resonance is considered to be in the range of 400-1200 MeV, while its width lies in 
the interval 600-1000 MeV [lo| . Clearly, one should not demand too much from the 't Hooft 
model with its simplifying assumptions. However, we consider it a significant finding that 
it predicts a low energy resonance at all, based solely on the leading-order quark exchange 
diagrams. 

In all four cases, the resonance is located very close to threshold, and it is very narrow. 
The width increases slowly with increasing quark-gluon coupling strength. On the other 
hand, the resonance position remains more or less unchanged as long as we restrict ourselves 
to models with real dressed masses (models I to III). Only for model IV, whose coupling 
constant is considerably larger, the resonance moves away from threshold to about 300 MeV. 

One might have expected the resonance energy to increase smoothly and in a more pro- 
nounced manner with the quark-gluon coupling strength g. However, one has to keep in mind 
that the included quark-exchange processes do not directly depend on g, but only indirectly 
through changes of the vertex functions and of the dressed quark masses that appear in the 
quark propagators. Owing to the already mentioned peculiar feature of the 't Hooft model 
that dressed quark masses decrease with increasing coupling strength g, the latter tend to 
decrease the resonance energy. Moreover, while a larger g implies an effectively stronger 
attraction between the two mesons - once they overlap - through the stronger quark-quark 
attraction, the very probability of this overlap drops in turn, because the spatial size of 



17 



the mesons decreases. These effects seems to counterbalance each other to a large extent, 
leaving the resonance position more or less unchanged. 

Similarly, that the width of the meson-meson resonance increases with g is also a conse- 
quence of the contraction of the mesons caused by the stronger quark- ant iquark attraction. 
This shrinking in size leads to a larger spreading of the meson-gg vertex function in momen- 
tum space (see Fig. [7J), which in turn contributes to the overlap integrals in the included 
Feynman diagrams in a wider momentum range, thereby broadening the resonance. 

Finally, we should mention that our principal finding - the existence of a narrow low- 
lying resonance in our calculations - does not depend on the particular value chosen as a 
constraint for the qq bound-state mass. If we use, instead of 140 MeV, a much larger or a 
much smaller value, we find again a narrow resonance close to threshold. 

V. CONCLUSIONS 

We calculated various models for quark-antiquark vertex functions within the 't Hooft 
model by solving the corresponding qq Bethe-Salpeter equation. In all cases, the bare quark 
masses and the quark-gluon coupling constants were tuned such that the mass of the qq 
bound state coincides with the pion mass. 

We found that, within a limited range of coupling constants and bare quark masses, one 
obtains qq bound states with the features of a constituent quark model, i. e., where the 
meson mass is approximately equal to the sum of the dressed quark masses. On the other 
hand, for larger values of the quark-gluon coupling constant this constituent quark picture 
is no longer sustained. 

We used the calculated Bethe-Salpeter wave functions to derive meson-meson ("pion- 
pion") scattering amplitudes within the 't Hooft model. They are calculated from the 
leading-order quark-exchange diagrams. These QCD-based meson-meson amplitudes exhibit 
a resonance structure close to threshold in all considered cases. We extracted an effective 
mass and width of this "sigma-like" resonance by comparison with a simple s-channel res- 
onance model where "pions" and "sigmas" instead of quarks are the effective degrees of 
freedom. 

The extracted values are not meant to be completely realistic in the sense that they 
should reproduce the experimental data, since the 't Hooft model is QCD only under sim- 



18 



plifying assumptions. However, the results of this calculation demonstrate that the 't Hooft 
model can accomodate a resonance for the four-quark system, already in the leading-order 
quark-exchange processes. In those exchange mechanisms, the diquark correlation, given 
by the quark-antiquark vertex function, plays an important role in determining the energy 
dependence of the meson-meson scattering amplitude. 

The fact that in our 't Hooft model calculations a narrow resonance lies close to threshold, 
while the broader sigma resonance is supposed to be located at higher energy, is indicative 
to the limitations of the 't Hooft model. It also hints at the importance of higher-order 
quark exchange processes as well as of gluon exchange contributions (higher order terms in 
the l/N c expansion), which should be investigated in future work. 
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